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1. Introduction 

The precise description of cooperative effects in strongly interacting many-body 
systems continues to pose many challenges. Paradigmatic examples are systems 
which may be described in terms of diffusion-limited reaction-diffusion processes. 
Applications of these systems and their non-equilibrium phase transitions have arisen 
in fields as different as solid-state physics, physical chemistry, physical and chemical 
ageing, cosmology, biology, financial markets or population evolution in social sciences. 
If the spatial dimension of these systems is low enough, that is d < d* where d* is the 
upper critical dimension, fluctuation effects dominate the long-time kinetics of these 
systems and their behaviour is different from the one expected from the solutions of 
(mean-field) reaction-diffusion equations, which attempt to describe the interactions 
of the elementary constituents in terms of the macroscopic law of mass-action, see e.g. 

[SZllinillllllMllIH]. 

One of the motivations for this work is the continuing practical interest in systems 
with reduced dimensionality, and such that homogenisation through stirring is not 
possible. Specifically, we shall consider the one-dimensional coagulation- diffusion 
process^ which is defined as follows. Consider a single species A of indistinguishable 
particles, such that each site of an infinitely long chain can either be empty or else 
be occupied by a single particle. The dynamics of the system is described in terms 
of a Markov process, where allowed two-site microscopic reactions A -I- o ^ -\- A 
and A-\-A^A-\-(ll or (ll-\-A are implemented as follows: at each microscopic time 
step, a randomly selected single particle hops to a nearest-neighbour site, with a rate 
r := Da^, where a is the lattice constant. If that site was empty, the particle is placed 
there. On the other hand, if the site was already occupied, one of the two particles is 
removed from the system with probability one. This model is one of the best-studied 
examples of a diffusion-limited process and at least since the work of Toussaint and 
Wilczek [46 it is known that the mean particle concentration c{t) ^ f^l'^ for large 
times and with an amplitude which is thought to be universal as confirmed by the field- 
theoretical renormalisation group j24[ [7]; in contrast a mean-field treatment would 
have predicted c(t) ^ t~^. These theoretically predicted fluctuation effects have been 
confirmed experimentally, for example using the kinetics of excitons on long chains of 
the polymer TMMC = (CH3)4N(MnCl3) [23], but also in other polymers confined to 
quasi-one-dimensional geometries [36l [21] , see also the reviews in [37] . Another recent 
application of diffusion-limited reactions concerns carbon nanotubes, for example the 
relaxation of photoexcitations [41 or the photoluminescence saturation [44|. On the 
other hand, the ID coagulation-diffusion process has also received attention from 
mathematicians [S] [21] and is simple enough that it can be related to integrable 
quantum chains, see [3l [42]. Hence, by a consideration of the quantum chain 
Hamiltonians, which can be derived from the master equation, the time-dependence 
of its observables could in principle be found via a Bethe ansatz, see [40]. In practice, 
however, it has turned out to be easier to find the time-dependent densities from 
the empty-interval method, which considers the time-dependent probabilities En(t) 
that n > 1 consecutive sites of the chain are empty [5j [TTl [4| |6| [27], see also [43] . 
The En{t) satisfy a closed set of differential-difference equations, subject to the 
boundary condition EQ{t) — 1 and the average particle concentration is obtained 
as c{t) = (l — Ei{t)) /a. The scaling behaviour of the averages can be directly studied 
in the continuum limit a — ^ 0, when En{t) — > E{x,t) which in turn satisfies the 
diffusion equation (9t — 2Ddl)E{x,t) = with the boundary condition E{0,t) = 1 
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such that the concentration now becomes c{t) = ~dxE{x,t)\^_^. Still, the direct 
solution of the problem is usually considered to be complicated enough to prefer to 
consider instead p{x, t) := dlE{x, t) where the boundary condition becomes p(0, t) = 
such that standard Green's functions of the diffusion equation can be used, see [5] and 
references therein. 

Remarkably, the empty-interval method can be applied to a large class of 
coagulation-diffusion models, where several additional reactions can be added, see 
e.g. [Sl[IIllllll271[llllI71[ini[ll3Il[31[2E]- Furthermore, the quantum 
hamiltonian/Liouvillian of coagulation-diffusion models with the (reversible) reaction 
2 A o A can, by a stochastic similarity transformation [22l [T5j [10] , be transformed to 
the one of pair annihilation/creation 2 A o which in one dimension can be solved 
by free-fermion methods, see e.g. [9l |45l [39l EH [3j [TH |43l [30]. Those one-dimensional 
reaction-diffusion systems which can be treated with free-fermion methods have been 
classified [16l |42], but the empty-interval technique has the advantage that further 
reactions can be treated, such as — > A or A(i}A — > AAA, which have no known 
analogue in a free-fermion description. In particular, Peschel et al. |3S] suggested a 
systematic way to identify observables for which closed systems of equations of motion 
can be derived from the reformulation of the master equation in terms of a Hamiltonian 
matrix in a controllable way. Their approach includes the method of empty intervals as 
the most simple special case. In principle, their method can be extended to include the 
probabilities of having several empty intervals of sizes m, n2, . . . at certain distances 
which allows to find correlation functions as well. Their study is the main subject of 
this paper. 

In particular, our approach allows to consider arbitrary initial configurations 
of particules and hence our results will include many of the existing results in the 
literature as special cases. As we shall see, there exists a natural decomposition of 
the time-dependent observables which may be arranged in terms of the information 
required on the initial state. This can be formulated through single-interval or two- 
interval probabilities for those quantities which we consider explicitly. We shall give 
examples which suggest a clear order to relevance in the long-time limit. On the other 
hand, we shall assume spatial translation-invariance from the outset, which simplifies 
the equations to be analysed. However, if one were to investigate the effects of disorder, 
one would have to revert to a formalism [111 |3] where translation-invariance is not 
required. 

The study of correlation functions of reaction-diffusion systems is also motivated 
by the recent interest in ageing phenomena: having begun in the study of slow 
relaxation in glassy systems brought out of equilibrium after a rapid change in the 
thermodynamic parameters, it was later realised that the three main characteristics 
of physical ageing, namely (i) slow, non-exponential relaxation, (ii) breaking of time- 
translation-invariance and (iii) dynamical scaling also occur in many-body systems 
which in contrast to glasses are neither disordered nor frustrated, see [l^ for a brief 
review and a forthcoming book 119] . Furthermore, these characteristics have also been 
found in several many-particle systems with absorbing stationary states, such as the 
contact process [12j |38l |8] , the non-equilibrium kinetic Ising model [33] or kinetically 
constrained systems such as the Frederikson- Andersen model [28]. One particular 
point of interest in these ageing systems is the relation between two-time correlations 
and responses and a study of the coagulation-diffusion process (along with its exactly 
solved extensions) should be useful, since exact results can be expected, at least in 
one dimension. However, while such an analysis is readily formulated in terms of 
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the empty-interval method at two different times, the explicit calculation requires 
the knowledge of the exact equal-time two-interval probabilities. In this paper, we 
shall provide this information, which will become an initial condition for the two-time 
correlator, and is going to be used in a sequel paper where the ageing behaviour in 
exactly solvable reaction-diffusion processes will be addressed. 

This paper is organised as follows. In order to make the presentation more self- 
contained, we recall in section 2 the derivation of the equation of motion for the 
empty single-interval probability En {t) before we proceed to show that the boundary 
condition Eoit) = 1 can be fixed through an analytic continuation to negative values 
of n. The techniques thereby developed are to be generalized to the two-interval 
probability in the remainder of the paper. The passage from the initial state towards 
the scaling long-time regime as a function of the initial distrbution is analysed and 
we also compare between the discrete model and its continuum limit. In section 3, 
the equations of motion and the formal solution is given, to be followed by the 
derivation of the consistency conditions with the single-interval probabilities. The 
general two-interval probability for arbitrary initial conditions is derived in section 4 
and in section 5 we use the results for the derivation of the equal-time correlators. 
We conclude in section 6. Several appendices (A-G) contain technical details of the 
calculations. 



2. Single-interval probability 

2.1. Equations of motion 

Using the definition of the coagulation-diffusion process as given in the introduction, 
we begin by recalling the derivation of the equation of the empty-interval probabilities 
[5]. The same equations can also be found within a quantum Hamiltonian formalism 
[35] . but this wiU not be repeated here. We denote by En{t) the time-dependent 
probability of having an interval of n consecutive empty sites at time t. Since the 
system is assumed to be homogeneous, En{t) is site- independent and will depend only 
on the interval size n and time t. The time evolution of this quantity is governed by 
the rate at which particles move on adjacent intervals of size n or n — 1. In an interval 
of length n, which will be denoted by | n | , a particle (•) can enter from the left or 
the right between the time period t and t + dt, and En{t) decreases during this period 
of time by the amount 



Pr(» ^ I n I) + Pr(| n | ^ •) 



The probability Pr(» ^ | n | ) is proportional to the probability that a particle lies on 
the left of the interval, or Pr(» ^ \ n | ) = Pr(« | n "| )rd/:, which can be evaluated using 
the relation 

Pr(>nn) + Pr(°E) = Pi-(E) = En{t) 



where the symbol (o) refers to an empty site. Since by definition Pr(o | n "| ) = En+i{t) 
we obtain directly 



Pr(»nn) = Pr(S.) = En{t) - En+i{t). (1) 

En{t) may also increase, if we consider the possibility that a particle sitting next to 
an interval of size n — 1 moves away from this interval. 



Pr(^ «|n-l|) + Pr(| n-1 |« ^) 
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This is possible because the process A + A ^ A constrains each site to contain at 
most one particle. Hence there is no need to consider the case when the particle 
encoun ters a nother p articl e when it moves away from the interval. As before, we have 
Pr(^ « |n-l I ) = Pr(« |n-1 \ )Tdt = {En-i{t) - En{t))Tdt. Summing the contributions, 
the rate of change for En{t) is given by 



dtEn{t) = 2T 



{En{t) - E^+i{t)} + {En-i{t) - E^{t)} 

^ 2r {En-i - 2E,, + E,,+i) . (2) 

This equation is valid only for a positive index n > 1. For n = 1 the rate of change 
for El (t) is given as previously by the equation 

dtEi{t)dt = Pr(» ^) + Pr(-^ •) - Pr(» o) - Pr(o ^ •) . 

We also have Pr(» ^) — Pr(»)rdi and Pr(» ^ o) = Pr(» o)rdt. The solutions for each 
of these quantities can be found by considering the probability conditions 

Pr(«) + Pr(o) = 1 ^ Pr(») = 1 - Ei{t) 

Pr(« o) + Pr(o o) = Pr(o) ^ Pr(» o) = Ei{t) - E2{t). (3) 
Therefore, the equation for n = 1 is given by 



dtEi{t) = 2T l-2Ei{t)+E2{t) 



(4) 



In order to be able to write this as the extension of eq. ([2]) for n = 1, it appears 
convenient and is, indeed, common, to introduce the constraint £^o(^) = 1- We shall 
do the same, but return to this condition below. However, the boundary conditions, 
including E{Q, t) — 1 were considered to be sufficiently complicated so that an explicit 
solution of ([S|) is usually avoided. Ingenious ways have been developed to extract 
physically interesting information, such as the particle-density c{t). We shall require 
the explicit form of E{x, t) below when looking for correlation functions and shall now 
give it. In the continuum limit, when a is small, we set x = na, and E{x,t) = En{t). 
The previous relation ([2|) can be expanded with respect to a and a rescaled hopping 
rate D = F/a^, which leads to a simple diffusion equation, together with a boundary 
condition 

dtE{x,t) ^2Dd,j,^E{x,t), and i;(0,i) = 1- (5) 

If we could use a spatially infinite Fourier transform E{x, t) ~ exp(ifcx)£'(fc, t) 

to solve the previous equation, we would obtain in the standard fashion 



E{x,t) = 



dx' 



exp 



— 7 (a; ■ 



x'f E{x',0), 



where the integrals over the real axis are unrestricted. In the above expression, a 
diffusion length 

io ■■= VSDt (6) 
acts as the scaling length of the function E{x,t) = E{x/£o). 



2.2. Effect of the boundary condition: continuum limit 

The simplistic approach outlined at the end of the previous subsection must evidently 
be modified in order to take the boundary condition E(0,t) = 1 into account. This 
amounts to define in eq. ^ the meaning of the probability E{x',0) for negative x' 
and is achieved by the following result. 
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Lemma 2.1 // one extends the validity of eq. (0) to all n G together with the 
boundary condition Egit) — 1, one has 

£;_„(i) =2-S„(<). (7) 

In the continuum limit, this leads to E(—x,t) = 2 — E(x,t). 

Proof: This is proven by induction. First, we consider the case n = 0. Using eq. ^ 
and Eo{t) = 1, we obtain 

dtEoit)^2T (E^i - 2Eo + El) , 

which impUes E^i{t) = 2Eo{t) — Ei{t) = 2 — Ei{t). In the general case, let us consider 
the equation of motion for the index — 1 and use the assumption ([7]) for the indices 
—n and — n + 1: 

E-n-l = 2E^n — E^n+1 + -^dfE^n 

= 2(2 - En) - (2 - Er,-i) + ^dt{2 - E„) 

— 2 — En+l 

where the equations of motion ([2]) were used again. This completes the proof, q.e.d. 

In the continuum hmit, this relation allows us to rewrite the integral ^ over the 
positive axis only 

/•oo A ^ r 

E{x,t) = erfc(x/£o) + / S(x',0) e-^"""')'/^"' - e-(-+-')V^o= (g) 

Jo VTT i:o L -I 



and where erfc is the complementary error function [T]. 

Eq. ([8|) is the general solution for the probability E{x, t) of having an empty 
interval, at least of length x and at time t, where the initial state is described by the 
function E{x,0). The particle concentration c{t) = Pr(»)/a can be obtained in the 
continuum limit from the relation ([3]): 

Pr(») + Pr(o) = 1 Pr(») = ac{t) = 1 - Ei{t), 

where ac{t) — 1 — Ei{t) ^ 1 — E{0, t) — adxE{x — 0, t), and therefore 

c{t) = -d,E{x,t)l^,. (9) 

The function E{x,t) can by definition be written as a cumulative sum of the 
probabilities for having bounded on the left, of size at least equal to x' or P{x',t) = 



Pr(. x' ) 



nOC 

E{x,t) = / dx' P{x',t). (10) 

J X 

This imposes two boundary conditions: first, we have E{0,t) — J^dxP{x,t) — 1 by 
normalisation. Then, in the limit x — >■ oo, one must have E(x, t) — > 0. 

We can express E{x,t) as function of P{x,0) by performing an integration by 
parts of dH): 

E(..t) . 1- l|"<i.'P(..,o)[orf(flt£) („) 

By differentiation with respect to x, we obtain the expression for the concentration 

c{t) = ^ / dx' P{x\ 0) exp ( - ) (12) 

From this, all initial conditions, characterised by P{x,0), lead to the long-time 
behaviour of the concentration: 
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Lemma 2.2 For sufficiently long times and any initial distribution P{x,0), the 
concentration decreases as 

c{t)c,^^+o{£o-') (13) 



Proof: Since P{x,0) is a normalised probability distribution, J^dx P{x,0) — 1, we 
must have P{x, 0) = o{l/x) for x oo. We rewrite eq. (|12p as follows 

9 1 9 /"OO 

c(0 = ^7 ^ dxPix£o,0)(l-e--') (14) 



o(l/fo) 

where the last estimate follows from the large- a; behaviour of P{x, 0). q.e.d. 

In a more explicit way, this may be obtained from eq. (I12p by a formal expansion of 
the exponential. If the second moment of P is well-defined, this leads to the long-time 
behaviour 

c(t)^^(l-^)-fo(£o-^). (15) 

On the other hand, the case of a diverging second moment is illustrated by the example 
P(x,0) = ^ °° , (16) 

where < a < 1 and oq = (1 + a) sin[7r/(l 4- a)]/7r is the normalisation factor. In this 
case, a calculation analogous to the proof of lemma 2.2 gives 



c(t) = / dxao- 



exp 



(-«i) 



2 



^4 



1 - ^-^ I dx ■ 



4" Jq 



(17) 



which gives the leading correction in the long-time limit as function of the 
exponent a. We also notice that if particles occupy each site with probability p, the 
concentration is defined by cq — p/a. When the system is filled with a concentration 
Co of particles, the function i?„(0) is proportional to 

En{0) ^ (1 - pT = (1 - ac^T''' E{x, 0) - exp(-cox) (18) 
However, if the system is entirely filled with particles, p — 1 and i?„(0) = for 
n ^ 0, then, from ([8]), E{x,t) is simply given by erfc(a;/^o), and c{t) = 2/y/TTio- In 
the general case of a given concentration cq where Eo{x) = e~^°^, we simply have 
P(x,0) = coe-'=o=^ and from 1^2^ 

1 2/, 2 \ _„r- ( 



c{t) = CO exp ^^4£o'j erfc j . (19) 

In figure [T] we illustrate the effect of several initial empty-interval distributions 
E{x,0). Clearly, as expected from the above discussion, all initial distributions lead 
to the same long-time asymptotics c(t) ~ but the way this asymptotic regime 

is reached depends on the initial state. This can be better understood when plotting 
the interparticle distribution function (IPDF) 

, , 1 d^E(x,t) 
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Figure 1. Time-evolution of the concentration c{t), for several initial conditions 
expressed in terms of E{x, 0) and parameter cq = 1/2. For any initial distribution, 
the particle concentration shows the same asymptotic behaviour. The passage 
between the initial and the asymptotic regime can be qualitatively explained in 
terms of the interparticle distribution function (IPDF) defined in eq. I I20II and 
shown for t = in the inset. See text for details. 

which gives the probabihty density that the next neighbour of a particle is at distance 
X at time t (SJ |B] . This function is shown for i = in the inset of figure [1] One 
observes that in those cases when p(x, 0) decays monotonously with x, the transition 
to the asymptotic regime is more gradual. On the other hand, in the third case there 
is an initial non- vanishing distance the particles must overcome before they can react. 
This leads to a very sharp transition between the initial and the asymptotic regimes. 

2.3. The discrete case 

We now give the solution of the discrete case, without performing a continuum limit. 

First we recall the solution of the differential equation The generating 

function F{z,t) = X^n^cx) '^"^n(^) satisfies as usual the differential equation 
dtF{z,t) = 2D(z + l/z-2)F{z,t), with the solution 

+ 00 +00 

n— — oo m— — OO 

Identifying En{t) in the previous expression, we write 

+ 00 

i^nW^e-*^* ^ £;„(0)/„_™(4Dt). (22) 

m — — oo 

We now must take the boundary condition EqII) = 1 into account. As in the 
continuum case, we replace the summation over negative values of the index m by 
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Figure 2. Time evolution of the concentration in the discrete case (full curves) 
and the continuous limit (dashed lines). The initial concentration cq is [1,0.5,0.1] 
from top to bottom. 



using the discrete relation £'_„(t) ~ 2 — En{t) and find 

"+CXD +00 

^ Ern{0)(ln~m{^Dt) - Ir,+m{^Dt)) + ^ 2/„+„(4i?t) + h{ADt) 
.m— 1 m — 1 

In the discrete case, the particle concentration is given by c{t) = 1 — Ei(t). Using 
summation and recurrence relations over modified Bessel functions, we obtain 

c{t) = [h{^Dt) + h{4Dt) - ^^™(0)^™(4^^) ) > (24) 

which generalises earlier results of Spouge [43]. In figure [21 we compare the particle 
concentration according to the discrete case eq. ([Mj) with the previously obtained 
solution eq. (jl9[) in the continuum limit, for three values of the initial concentration 
Cq. We used, respectively, the initial distributions £'„(0) = (1 — co/a)" and E{x,0) ~ 
g-cox^ As expected, the same asymptotics is found in all cases, independently of 
the initial concentration. We observe that the passage between the initial and the 
asymptotic regimes is more gradual in the continuum limit. As above, we interpret 
this as coming from the fact that in the discrete case the particles must first overcome 
a finite distance before then can react. 



En{t) = e~ 



3. Tv^ro-interval probability 



In this section, we generalize the previous result and evaluate the probability 



Eni,n2{d,t) to have two empty intervals, at least of sizes ni and n2 and separated 

Prf 



by the distance d: we denote it by E^ ^2 ("^i ^) 
expected to have the following symmetries 



ni 



"2 



This function is 
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Eni ,0 (d; t) ^ (t) and £'o,„2 {d, t) = E^^ (t) 

Eni,n2{0,t) = i?„j+„2(<). 



(25) 



3.1. Equations of motion 



As before and using the notations of previous section, we consider the different 
possibihties for the variation of E'n^ ^^C'^:^) between the time t and t + dt: 



+ Pr( 



ni 



n2 



ni 



d 

^ •d- I 



Pr( 



ni 



n2 



n2 



Pr(^ • ni — 1 
- Pr( 



ni — 1 



d 

• ^ d 



n-2 



+ Pr( 
-Pr 



ni 



d — 1» ^ 



(26) 



n2 



ni 



712 — 1 • ^) 



n2 



)+Pr( 



ni 



d ^ 



?l9 — 1 



The probabihty rates are given by considering the sum rules for static probabihties. 
First, we consider the negative contributions for which we obtain the relations 



ni 



"2 



= Pr(»| ni d 



Pr(« 



Pr(»| ni I d I n2 ) + Pr(o| ni d 



n2 



^^2 



)rd< 

Pr( 



ni 



^^2 



Pr(»| ni d 



n2 



) — £'ni,n2 (d, t) — Eni+l,n2 (d, t) 



and 



Pr( 
Pr( 



ni 



ni 



^ •d—l 
• d-l^ 



"2 



n2 



) = Pr( 
HPr ^ 



ni 



id- 1 



"2 



Ofi- 1 



n2 



)rdi, 
Pr(" 



n2 



Pr( 



id- 1 



n2 



For the positive contibutions, we have 
Pr(^ • ni - 1 rf I rt2 1) = Pr(< 



rti — 1 



)rdi, 



"2 



Pr(» ni — 1 



«2 



Pr(o ni — 1 



n2 



Pr( ni — 1 



f^2 



Pr(. 



ni -=4 



^^2 



) — £'ni-l,n2 (c^j ^) ~ -^ni,n2 (^^j 0: 



and 



Pr( m - 1 



Pr( ni — 1 



• ^ d 

d^ 



«2 



Pr( ni - 1 



• d 



n2 



Fdi, 



n2 



) + Pr( ni - 1 



o d 



f^2 



Pr( ni — 1 



«2 



Pr( 



ni — 1 



• d 



n2 



) — Eni-l,n2 {d + 1, i) — Eni,n2 {d, t) , 



(similarly for the other terms which are symmetric). After gathering all the 
contributions, we finally find (the time variable is from now on suppressed) 

ati;„i,„2W =r[-8i;„,,„,(d) (27) 

+ En 

1 +1 ,n2 

{d)+E 

ni,n2 + l(^) + ^ni-l,n2(^) + ^ni.n2-l(^) 
(d - 1) + Eni^n2 + l{d - 1) + £'ni-l,n2(^ + 1) + Eni,n2-l{d + 1)] • 

We have checked that the same closed system of equations of motion is also obtained 
when the master equation is rewritten in terms of a quantum Hamiltonian [35] . 
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The continuum limit of this diffusion equation is obtained by expanding the terms 
up to the second order in the lattice step a (the will be absorbed in F). Setting 
X = nia, y = and z — da we obtain the following linear differential equation: 



dtE{x,y,z) = 2D 



dl + dl + dl 



dxdz + dydz 



E{x,y,z) 



(28) 



3.2. General solution 



The general solution for (^5]) is obtained by diagonalising the quadratic form associated 
with the differential operator P := d'^-\-dy-\-dl— dxdz— dydz. We find that the following 
change of variables (x, y, z) (X, Y, Z) diagonalises P: 

X = a{x + y^^z) , y = /3(x + y - v^z) , Z = 7(2; - y), (29) 

with the positive constants 

a = [2(2 - \/2)]-i/2 ^ /3= [2(2 + V2)]-i/2 , 7=1/^2, 

such that + = 1 and - 0^ ^ 1/V^. 

In the new variables, the operator P — d\ ^ d\ ^ d\ \s diagonal. As in the 
previous section (see ([6])) for the one- interval problem, and in the continuum limit, if 
the boundary conditions are ignored for a moment, the function E{x, y, z) can be found 
via a Fourier transformation, and explicitly expressed as a kernel integral depending 
on the initial conditions Eq{x, y, z) := E{x, y, z, 0) : 



E{x,y,z) 



V2 



dx'dy'dz' W{x-x',y-y',z^ z')Eo{x', y' , z') (30) 



where the Gaussian kernel 'W{u,v,'w) is given by 

yV{u,v,w) = exp^ - [u + V + \/2wY - {u + V ~ V2w)'^ ~ ^'^ [u - v) 
io L 

1 



1 r 

to 



{u + V + w) 



2 2 
w 



-{u - vY 



(31) 



the Jacobian of the transformation ((29)) being equal to 4:\/2af3j — \f2. 



3.3. Compatibility conditions 

As expressed before by equation ([7]), it is important to make a correspondence 
between intervals of formally negative and positive lengths. This will be needed 
in the formal solution (|30|) which requires real variables, whereas the probability 
E{x, y, z) has an obvious physical meaning depends for positive distances only. We 
have to consider 3 cases, depending on whether x, y oi z are negative or positive. By 
symmetry considerations p5|). it is only necessary to consider the case where x or z 
are individually negative (the case y < being deduced by means of the first equation 
([25])), the other variables being positive. The explicit evaluation of eq. ((30|) in the 
following sections requires several identities, stated as lemmata for clarity and proven 
in appendix A and B, respectively. The first one treats the case of formally negative 
interval lengths. 
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Lemma 3.1 The probability of two- empty-intervals of negative lengths is related to 
the probability of positive lengths as follows. 

Eni.-riiid) =2Eni — i?rii ,"2 (^^ ~ "2 ) , (32) 

E^n^.-n2id)=4-2En,-2En2+En,,n2id-ni~n2). (33) 

A further relation connects the negative separations — d between two intervals to the 
positive ones. 

Lemma 3.2 We have 

-E„j,„2( — d) = 2Eni+n2-d — En^^d,n2-d{d)- (34) 

Later on, we shall require these results in the continuum limit, where the 
expressions ([7]), ((32)) . (|33l) and (l34l) take the following form 

E{-x) ^2-E{x), 

E{-x, y, z) = 2E{y) - E{x, y,z-x), 

E{x, -y, z) = 2E{x) - E{x, y,z- y), 

E{-x, -y, z)=4- 2E{x) - 2E{y) + E{x, y,z-x-y), 

E{x, y, —z) — 2E{x + y — z) — E{x — z,y — z, z). (35) 

This allows us to rewrite pop in the restricted domain where {x' , y' , z') are all positive, 
and where Eq{x' , y' , z') is physically well-defined. 

4. General solution for E{x,y, z,t) 

From the general equation ([50)1 . we separate the 8 different domains of integration 
around the origin for example (x' > 0,y' > 0,z' > 0), {x' < 0,y' > 0,z' > 0) etc., 
and use relations (O, ([5^ . ([55)) and (|M)) to map all domains into the single domain 
(x' > 0,y' > 0, z' > 0). This calculation is done in the appendix C and, here, we just 
summarize the results. The general solution can be decomposed as follows: 

E{x, y, z, t) = -B^o) (x, y, z, t) + E^^'^ (x, y, z, t) + E^^^ (x, y, z, t), (36) 

where E^^"^ {x,y, z,t) is obtained from the terms independent of the initial conditions, 
E^^"^ (x, y, z, t) from the initial one-interval probability Eo{x') and i?*-^-* (x, y, z, t) from 
the initial two-interval probability Eo{x,y, z), respectively. Note that E'^^\x,y, z,t) 
and E^'^\x,y,z ,t) depend on Eq{x'), EQ{x',y',z') with arguments positive, hence 
this gives us the physical answer to the diffusion process in the coagulation problem 
starting from arbitrary initial conditions and constraints on the differential equation. 
We now analyse these three terms one by one. 

Special case of a system initially entirely filled with particles 

We notice that (IC6p . ()C7p . and ()C8[) contain initial conditions for the single- interval 
distribution Eq{x'), some constants independent of the initial conditions, and initial 
conditions for the two-interval distribution Eo{x' ,y' , z'). To simplify notations, we 
shall re-scale aU lengths by Iq such that E£g{x,y,z) = E{xio,y£o, z£o,t). In ()C8p . we 
can isolate from Ea{—x', —y', z') two terms independent of the initial conditions, 

4 - Ae{y' - z')e{x' - z')e{x' + y' - z'). 
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It is obvious that e{y' ~ z')9{x' - z')9{x' + y' - z') = e{y' - z')e{x' - z'). These two 
terms, plus the first term in (jC3l) . give a contribution to the general function equal to 

Ef'J (x, y,z) = crfc(z)erfc(a; -\- y + z) 

+ [ dx'dy'dz' m^i-x',-y',z'){4- 49{y' - z')eix' - z')}. 



■ JrI 

Performing the translations x' — z' x' and y' ^ z' ^ y' in the last contribution, we 
obtain 

E^^^ (x, y, z) = erfc(z)erfc(a; + y + z) 

+ J\ [ dx'dy'dz' \m„{^x',-y',z')-m„i~x' ~z',-y' -z',z')}. 

Using the relation (jC4p . and the identity 

(x^x' -z',y-y' ~z',z + z') = We„ {x ~ x' ,y ~ y' , z ~ z')e~^''' 
^W,,{x' +z\y' + z\-z') = -W,,{x',y\z') (37) 
we can rewrite E as 

Ef^ {x,y,z) — erfc(z)erfc(a; + y + z) 



+ \^ da; 



+ 



,'dy'dz' {Wioi-x', -y', z') + W,o(-x', -y', -z')} 



erfc(z)erfc(a; + y + 2) 



dx'dy' / dz'Wt,{-x\-y',z') 



The integral over z' now gives a gaussian exponential. Then the two remaining 
integrals can also be carried out explicitly (see appendix G). Introducing again the 
diffusion length £0, we find 

S<»)(.,.,.) . »tc (I) crfc (1) + „r,c (i) (i±|±i 

_erfc(i±i)o,fcf=^). (38) 



to 

This is the exact two-interval probability in the case of an initially fully filled lattice 
(where both Eo{x) and Eo{x,y,z) vanish). In particular, the solution psp satisfies 
the symmetry conditions (j25l) . In the limit of z large, one has the factorisation 
Eix,y,z,t)^ E{x,t)E{y,t). 

The remaining terms of the full solution depend on the initial conditions. They 
are of two kinds, and involve either the single-interval or else the two-interval initial 
probabilities. We turn to them now. 



4-2. Contributions to E[x^y,z,t) from terms with a single- interval initial distribution 

The contributions to E{x, y, z, t) of single- interval distributions come from the previous 
relations (jC6|) . (jC7|) and (jC8[) . where we can isolate the following individual terms 

• the second term of equation (|C3|) 

• the first 3 terms in ([06]) and JCTl) 
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On the whole, there are 10 terms contributing to the initial conditions given for a 
given choice of Eq[x'). Gathering these terms and performing successive translations 
in a;', y' or z' when necessary, we obtain 

«(x,y,z) = erfc(z) T ^ E^^,,{x')(e-^-+y+'-'^"^' - e-(W-+-')^ 



+ / dx'dy'dz' Eo{x')\m„{x',~y',z')~m„{^'^z',-y'^z',z') 

V TT'^ Jr3. ' 

+ m„ {-y\ x', z') - i-y' - z\ x'~z', z') 



+ W,„(-x' - z', -y' - z', z') - m„ i^x', ~y', z') 



+ m„{-y' - z', -x' - z', z') - \Vtoi-y', 



+ e{-x' -y' + z') [W,„ [y', ^z', x') + {~z' , y' , x') 



+ e{~x' + y' + z')W,„{-z',-y',x')] 



(39) 



By performing partial translations and integrations and simplifying all terms from Ji 
to /g (see appendix D for details), we obtain 

dx' „ ^^-(x'-x-y-zf _ ^-{x'+x+y+zf 



Eel\x,y,z)= / Eo^^„{x')l^ciic{z) 



+erfc(a;) 
+erfc(y) 



^-{x'-yf _ g-(an'+?,) = 



^-(x'-xf _ ^-(x'+x)"' 



erfc(a; + z) 
erfc(y + z) 



^-{x'-y-z) _ ^-{x'+y+zY 



^~(x' -x-zY _ ^-(x'+x+zf 



+erfc(a; + y + z) 



^-(x'-zf _ ^-(x'+zf 



°° Ax' 

—r= Eo,io{x')Ki^ia{x';x,y,z), 



(40) 

where the kernel Ki,ig is positive. When z = 0, we recover the result (|8]) for a single- 
interval distribution of size x -\- y. For some functions Eq(x'), the previous integrals 
can be performed exactly since Ki is gaussian in the variable x' . For example, if we 
take as initial function Eq{x) = e"*^"^, where cq is an initial concentration of particles, 
we find 

y, z) = erfc F,„{x + y + z) + erfc (^|) F,,{y) - erfc (^^) F,„{y + ^) 



erfc ( — I Fcg (x) — erfc 



y + z 
4 



Fco(x + z) + erfc 



X + y + z 

To 



with the following abbreviation 
dx' 







^-(x'-xf/t^ _ ^-(x'+xf/l,,^ 



Fcoiz), (41) 



(42) 
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The limiting values of this function read 
and Fq{x) = erf(a;/£o)- 



In the long time limit, the above expression goes to zero like l/^o 



Ax 



1 



6 



cg4' lo 



(44) 



and the dominant part of contribution E^^\x,y, z,t) in the same limit behaves like 

(45) 



It is interesting to compare this expansion with the long-time limit expansion of 
E^^^^x^y, z), which is independent of cq 



(46) 



The first term (|45|) tends to increase the two-interval probability by a factor 
independent of the distance, since there are less particles in the system for a finite 
concentration of particles. 

^.3. Contributions to E(x,y,z) oj two-interval initial distributions 



As noticed previously, we can isolate from equations ()C6p . \C1\ . and (jCSP terms 
involving E{x' , y\ z', t ^ 0) = Eo{x' , y' , z'). In particular 

• 1 term EQ{x',y',z') when all variables are positive, in combination with 

• 2x3 terms in JCU and JCll) 

• 5 terms in ((C8)) . 

When combining and simplifying these terms (see appendix E for details), we obtain 

(2) 

a reduced integral form for Ej, {x, y, z), as function of a kernel K2^io{x' , y' , z'; x, y, z) 



^fo V^^)^ \ ~3 f da;'dj/'dz' Eq/o {x' , y' , z')W^o - ^' ,y ~ v' , ^ - z')K2,£o {x' , y' , z'; x, y, z) 



with 



K2j„{x', y', z'; X, y, z) = [1 - e-4(.'+y+.')(W^)](i _ ^-^^'^^ 

^■-iy'y^l _ g-4z'(a;+j^+z)jj-[^ _ ^-4{x' +y' +z')zj 



-4x 



^-i{y'+z'){x+y+z)ll^ _ -A(x'+z')z 



-4(x'+z')(x+y+z)^ 



x—Az' [x+z) r 



-iy'y-iz' iy+z) r 



^^-iy' (x+y+z) 
^—Ax' (x+y+z) 



](1 
](1 



-i(y'+z')zT^ 
-ix' z\ 



(47) 



It is important to notice for the following sections that _ft'2(a;', y', z', 0, 0, z) 
dxK2{x' ,y' , z' ,x,y, z)\^^y^^ = dyK2{x' ,y' , z' ,x,y, z)\^^y^^ = 0. 
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4.4- Sum rules and large- distance limit 

From the resuhs of the previous subsection, we can put the total two-interval 
distribution as a sum of 3 terms 

Eix,y,z,t)^E^°^ {x,y,z,t) (48) 
dx' 



V2 



£o'7r3/2 



Eo{x')Ki{x';x,y,z) 
dx'dy'dz' Eo{x\ y', z')Wix -x',y - y', z - z')K2{x' ,y' , z'\ x, y, z). 



+ 

As a check on our calculations, we consider the particular case of initial conditions 
given by a configuration where there is no particle in the system (or cq = 0). In this 
case, E{x) and E{x,y, z) are always equal to unity for any time. This means that if 
we put the conditions Eo{x') = 1 and Eo{x' ,y' , z') — 1 into equation we should 
recover the result E{x,y,z,t) — 1. The first contribution comes from ([55]) and is 
independent of initial conditions. The second contribution comes from (|40p or (|4T|) 
with Cq — 0: 

f x \ I y \ /"^ y ' 

E{x, y, z, t) — erfc I — 1 erfc I I + erfc — I erfc 



0/ V^o/ V^O/ V t^O 



X + 




erfc 








erf 




X + 




erf 


4 





-erfc ( — - — I erfc ( ^ ^ " I -|- erfc I — I erf (4-1-1- erf I — 1 erfc i 4- 



4 J \ioJ VoJ \ioJ Vo 

-|-erfc [ — I erf [ — — r~- — I + erf ( — ) erfc J ^ 



io J VqJ V io 
.eric ( ^ I erf (' - erf ( ^] erfc ( ^+ ' 



V2dx'dy'dz' 



io J \ io J \ i 



W{x -x',y- y', z - z)Ki{x' , y\ z' \ x, y, z). 



The terms can be rearranged and we obtain the following equality, since E{x, y, z, t) 
1 for all times, 

V2 



v 2 f 

-3—— / dx'dy'dz' W{x-x',y-y',z-z')K2ix',y',z';x,y,z) 
TT'^/'^ Jr\ 

: erf ( ^] erf ( + erf ( ^] erf ( ^llll] erf ( ^] erf ^ ^ 



(49) 



JoJ \ioJ V i 

The last expression gives an identity for the complicated integral involving only the 
different weights when the initial functions are set to unity. When compared to 
expression we see that this is the same expression except for the erfc functions 
are replaced by erf. We also see that all correlators vanish on the empty-lattice, 
as expected. From the general expression (|48|) we can compute the limit when the 
distance z is large. As seen previously with (p8|) when the system is entirely filled with 
particles, the quantity E^^\x,y,z,t) is easily factorised as lim.z^i E^^\x,y, z,t) = 
erfc(a:/fo)erfc(y/^o) — E{x,t)E{y,t). In the general case, for a finite concentration at 
initial time for example, we can show that the limit is still valid from (l48t . Indeed 
the two other parts E^^^^ and E^^^^ have simpler behavior in this limit. For E^^\ the 
kernel Ki given in (j40p has two dominant terms 

-{x'-x}^ _ ~{x'+x) 



Ki^i^{x';x,y,z)c:^eiic{x) e'^^'-^) _ e-("='+^) + ciic{y) 



e 
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E'^^\x,y,z,t) ~ erfc(x/4) 



dx' 



E{x',0) 



dx' 



+crfc(y/4) 
The other kernel K2 behaves hke 

K2,ioix',y',z';x,y,z) ~ 



E{x',0) 



(50) 



1 - e 



1 -e 



-4yy' 



Then 



E<-^\x,y,z,t) 



V2dx'dy'dz' 



Eoix ,y ,z )yViaix ~ x',y ~y',z - z') 



x(l~e-''^^')(l-e-'^™') 



In the last expression, we set z" = z' — z, so that Eo{x' , y' , z') — Eq{x' , y', z" + z). We 
also assume that Eo{x' ,y' , z" + z) ~ EQ{x')Eo{y'), and that the integral over z" can 
be extended over the real axis since the lower bound —z is large and negative. Then 
it remains an integral over the weight Wio {x — x' ,y — y' , ~z"\ which is independent of 
z, and which can be performed exactly. We find in particular that 

^/2dz" 



yV,,{x -x',y~ y', -z") = e~{x-xr/i.'-{y-yr/^.\ 



If we multiply this exponential with the kernel K2, we obtain that 

dx' 



E^^\x,y,z,t) 







E(x',0) 



(51) 



dy' 
lF£c 



Eiy',0) 



-(y'+y)Vfo' 



The sum of the expressions (|50)) . (|51l) . combined with the limit of ((38)) . gives 
exactly, after factorization, the expected limit, which is the generalization of the result 
seen with ([55)1 in the particular case of an initially filled system 

XvaiE{x,y,z,t)^E(x,t)E{y,t) (52) 

2:>1 

This limit is exact for any given initial condition at any time. 

Summarising the contents of this section, we have the decompositon p6p of 
the two-interval probability, where the individual terms are given by eqs. p8l40l47p . 
respectively. 



5. Two-point correlation function 

5.1. Definition 

The two-point connected correlation function can be defined in the discrete case as 
the probability to have two particles separated by the distance d 

Cdit) := Pr(. d .) - Pr(.)Pr(.). (53) 

We can express the previous function in terms of the discrete two-interval functions. 
Indeed 

Pr(« d»)+ Pr(o d •) = Pr(«) = 1 - Pr(o) = 1 - Eo,iid). 
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Since Pr(o d •) + Pr(o d o) ^ Pr(o) = Ei^oid), we have Pr(o d •) = Ei^oid) - 
and the correlator becomes 

Cd{t) = 1 - EoAid) - E^^oid) + Ei,i{d) - (1 - Ei,o{d)f. 

Using the fact that Eo,o{d) = 1, we can expand Cd{t) up to the second order in the 
lattice step a, setting z — da, and C{z,t) — Cd{t)/a^, so that 

C{z,t) = d%E{x,y,z)\^^^^^^^~ d.M^)\^=^dyE{y)\^^^. (54) 

This is the general expression of the correlation functions, which depends on the one 
and two-interval probability functions. 



5.2. Decomposition : general formalism 

There are three contributions to C{z,t) which come from the second derivatives 
with respect to x and y of E'^^^ and i?^^^ respectively. From the two-interval 

probability eq. p6l) . we have the decompostion 

C(z,t) = C(°)(z,t)-HC(i)(z,t) + C(2)(z,t), (55) 
where the different contributions are 

C(°)(z,t) =52j,£;(o)(x,2/,z,i)U=o,y=o- (a,erfc(x/4)U=o)' (56) 

(-2, i) = dlyE^^^ {x,y,z, t) U=o.y=o - 2dJ{x)\^=o dyeiiciy / eo)\y^o 
C(^) {z,t)^ dlyE^^^ {x,y,z, t) U=o.,=o - {dJ{x%^o?, 



with 



/to' 



._,-(x' +xf I ta- 



rn 



From ([55]). we obtain the first non-connected contribution of C'-'^^(z,t), which 
does not depend on the initial conditions 



dlyE'^''\x,y,z) 



2sinh(zV42) + V^— erfc 
to 



and so 



C(°)(z,t) = 



4e- 



V^^erfc ( ^ 



7^/0(2/4) 



(58) 



/o(3) 



4 

TT 



(59) 



where the scaling function /o has the following limit behaviour 

; for 3 ^ 
; for 3 ^ 00 

Eq. ([58)) is in exact agreement with the result announced in the literature [Sj H] . The 
asymptotic forms ([5^ have also been obtained several times, either for the coagulation- 
diffusion process [28l |31] |32] or for the equivalent [22l [15] pair-annihilation-diffusion 
process, see [3l [TH [30l [42] and references therein. However, the present discussion is 
not restricted to the rather special case of an initially fully occupied lattice. As we shall 
see, the second and third contributions in ([55]) are corrections to the leading behaviour 
in the long time limit (see appendix F for the details of the proof). We obtain a 
hierarchy in the inverse powers l/£o, where |C^°^(z,i)| \C^^^{z,t)\ |C'^^'(2;,t)|, 
with a dominant contribution of order Xjl^^ , l/^o^ and l/i!o^, respectively. 
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5.3. Application 

In what follows, we consider the special case of initial conditions Eq(x') — exp(— cqx') 
and Eq{x' ,y' , z') = exp(—co{x' + y')), where cq is the concentration of uncorrelated 
particles at time t = 0. We assume that initially the two-interval distribution 
is independent of the distance between the intervals. It satisfies nevertheless the 
condition Eq{x' ,y' , z' = 0) = Eq{x' +y'). The solutions for the correlation function 
given in appendix by (|F1[) and (jF2p can be computed except for the triple integral, 
where the successive gaussian integrals cannot be expressed, at least to our knowledge, 
in terms of known special functions. Since we are merely interested in the long-time 
limit and in the influence of the initial conditions in this limit, we can nevertheless 
perform an expansion in x' /£o and y' /io inside the kernel K2,eo derivatives and the 
weight function (— x', z — 2') in (jF2p . This is so because only small values 
of x' and y' are relevant when cq^o is large when these integrands are combined 
with the exponential factor exp(— co(a;' -I- y')). The latter function renders the 
integral finite after expansion of the weight and kernel. Indeed the natural small 
parameter of the expansion is 1/ (co^o)i relatively to the other dimensionless parameter 
z/£q, and the series expansion is assumed to break down only in the limit of low 
concentration or short times. But the previous result (IF2[) provides a general form 
suitable for series expansion of the correlated function with generic initial distributions. 
The other connected part of the correlation function C'-^-'(z,i), which is defined in 
equation (|Fip , can be computed exactly with the chosen exponential initial condition 
Eq[x') = exp(— cqx'), or by deriving twice with respect to x and y the previous result 
given in (|4ip . In the asymptotic limit anfn for cqIq large, we obtain 



C(i)(z,t) = ^|/^2 f-3 + 24)erf(z/4)e-^'/'°' 



1 - i^V'"" + ^ . 15 - 20^ +4^ erf(./4)e 



(-3 + 9i^-2i^)e--V.oV,(l/(,^,„)4), 



We notice that these first terms contribute to the correlation function in the large 
time limit at least like l/€o^ times a scaling function of z/Iq. This is a correction to 
C'°''(z,i) which behaves like l/^o^ instead. 

The last term C(2)(z,t) can also be expanded as a power series of l/(co^o) whose 
coefficients are scaling functions of z jl^ 

C^^\z^t) = ^exp(-2.VV)f:%^ (61) 



where the first few polynomials P2k{z) read 

P^{z) = ~l~2z^ , Peiz) = 12 + 24:z^ ~ 16z^, 
Ps{z) = -156 - 312^2 -I- b28z^ ~ 96z^ 

Pio(z) = 2400 + 4800^2 - 15360z'' + 58882^ - 512^^ (62) 

In general, C{z,t) can be expanded as a series in l/(co^o)i with coefficients being 
scaling functions of z/£o. Only the contribution C^°^(z,t) is dominant at zero*'' order 
in l/(co€o), while C^^\z,t) contributes at order l/(co^o)^- Otherwise the dominant 
part of C(^)(z, t), as seen in equation ([6T|) is of order l/(co^o)^, therefore smaller than 
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\ 




Figure 3. Time evolution of the connected two-point correlation function C{z,t) 
for CO = 1 and z = 1/2. The full solid line shows the leading contribution 
C(0)(l/2,t), the dashed line also includes the effect of the leading correction 
C(°)(l/2,t) -|-C(l)(l/2, t) while the dashed-dotted line includes all contributions 
C(0) (1/2, t) + CW (1/2, t) + C(2) (1/2, £). 



the previous ones. The correlations are negative over the whole range of considered 
values of z. 

In figure[3l the time-evolution of the correlator C(l/2, t) is illustrated, for initially 
uncorrelated particles with a concentration cq = 1 and 8D = 1, so that £q — \/t. In 
particular, we plot, respectively, the leading contribution for i — > oo, C^*'^ (1/2, i), 
along with the curve resulting when the first correction C'^^-'(l/2,t) is included, as 
well as when both corrective terms C*-"'^-'(l/2, t) + (7*^^^(1/2, t) are added. In agreement 
with our asymptotic estimates, we clearly see that even for relatively small times, a 
clear hierarchy |C(°ni/2,i)| » |C(i)(l/2, <)| » |C(2)(l/2,t)| emerges. 

5.^. Extensions 

Having found the single- and two-hole probabilities, we can immediately derive further 
quantities of physical interest. For example, the effective reaction rate is controlled 
by the pair probability Pr(»») of finding two particles on neighbouring sites. Similarly, 
one may define the triplett probability Pr(» • •). Simple enumeration leads to 

( \ p OP , p ''^o 29^E{x,t) 

Pr(»») = E0-2E1+E2 = a 



Pr(...)= E0-3E1 + 2E2- E3 + Ei^i{l) 



x=0 



a-.0 ^3 



d^y^E{x,y,z)\^^y^^^^- d^^^E{x,y,z)\^y^^^Q (63) 

While the pair probability can be expressed in terms of the single-hole probability 
alonelil the triplett probability already depends on a two-hole probability as well. 

II The x-dependence of d^E{x,t)/c{t) describes the interparticle distribution function and has been 
analysed in detail in the past, see [S] [6] and references therein. 
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In the continuum limit, the pair probabihty reads, for the three examples of initial 
distributions already considered in figure [T] 



Pr(»») = < 



co^c 



Co 4 



Co exp erfc — - ; if Eoix) 



dx e 



(1 + coioxY 



; if £;o(a:) = (l+cox)-^ 



(64) 







if Eq (x) — erfc 



Co a; 



where Cq characterises the width in the initial state. In the first two cases, the long-time 
behaviour is given by Pr(..) ~ (2/VF)(co/4) i"^/^ when co is kept fixed. Although 
the long-time behaviour is algebraic, the associated amplitude depends explicitly on 
the initial distribution, in contrast to what we had seen in eq. ()13p for the particle 
concentration c(t). 

Similarly, from the results derived in this paper further correlators can be directly 
obtained, for example 



Pr(» • d» 

Pr(»»d»»)= Pr( 



= Eo~ 3Ei +E2+ -Bia(d) + Ei,i{d + 1) - E2,iid) 
dxyzE{x,y,z)\^ 

(X,y, ^)\x,y=o\ 

• •d»)-Ei+E2+ Ei,i{d + 1) + Ei^iid + 2) 
- E2,i{d + 1) - Ei^2{d + 1) - Ei,2id) + E2,2{d) 

dxxyyE{x, y, 2) Ij. j^—o ~ dxxyzE{x, y, 

dxyyzE{x, y, z)\^ y^g dxyzzE{x, y, z)]^^^^^ 



a->0 4 

= a 



(65) 



6. Conclusions 



We have analysed a method for the computation of time-dependent correlators in 
the ID coagulation-diffusion process. The relevant quantities are the probabilities of 
finding either a single empty interval of a given size or else two empty intervals of given 
sizes and at a given distance from each other. These probabilities satisfy simple linear 
differential or difference equations, but working out the explicit solution is rendered 
difficult through boundary conditions, whose symmetry properties are not the same as 
those of the differential equations. Rather than circumventing this difficulty, we have 
shown how, by analytic continuation to negative sizes and distances, the full solution 
may be found for an arbitrary initial distribution which in turn can be characterised 
in terms of single- and double-interval probabilities. 
Specifically, we have seen the following; 

(i) The leading long-time behaviour of the particle density is explicitly confirmed 
to be independent of the initial conditions, as expected from the field-theoretical 
renormalisation group |24| [7] and in agreement with earlier calculations using 
a formulation of the empty-interval method where spatial translation-invariance 
is not immediately built in [5j [6] Our results also allow to compare directly 
the continuum limit with the result found for the discrete lattice and we can 
derive systematically the finite-time corrections to the leading dynamical scaling 
behaviour. 

^ As discussed in section 5, the various asymptotic estimates of either the diffusion-coagulation 
process or the equivalent pair-annihilation-diffusion process are fully reproduced. 



Coagulation- diffusion process in ID by the empty-interval method 



22 



(ii) The leading long-time behaviour eq. (j58p of the density-density correlator not 
only agrees with the earlier asymptotic estimate [1] but is furthermore explicitly 
shown to be independent of the initial conditions as well. The corrections to the 
leading scaling behaviour have also been analysed and the relative importance of 
the initial single- or double-interval probabilities can be quantitatively studied. 

(iii) By considering general initial conditions, we have identified a natural 
decomposition, see eqs. (I36l55p . into terms which develop a clear hierarchy in 
the long-time limit. 

This information will become important in a sequel article where we plan to analyse 
the ageing behaviour in the coagulation-diffusion process or generalisations thereof, 
where methods analogous to the ones developped here should apply and the single- 
time correlators studied here will serve as initial values for the two-time correlators 
whose dynamical scaling behaviour will be searched for. 
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Appendix A. Proof of lemma 3.1 

Starting from the discrete case (j27p . we consider the time differential equation where 
ni = 0, which leads to 

dtEo^nAd)^dtEr^,. (Al) 

The first term is evaluated from (P7|) and is equal, after some algebra, to 

dtEo,n2{d^) — r 2i?„2_i — 4:En.2 + 2En2+l 

id) + (rf + 1) - 4S„, + Ei^n2 (d) + Ei,n2 {d - 1) 

Using ([2]), the right-hand side in (|Aip is equal to 

dfEn^ — 2T En2-i — 2£'„2 + En2+l ■ 

By comparing the last two expressions, we obtain the equality between ui — —\ and 
ni ~ 1 

[d) + E.i^n2 {d+l)+ Ei^n2 {d) + £^i,„2 (rf - 1) = ^En2 ■ 

We then differentiate this relation with respect to time, and use (P7)) for each of the 
4 terms on the l.h.s (also ([2]) for the evaluation of dtEn2) in order to obtain a new 
relation between index ni = —2 and index ni = 2 

E^2.^2 (d) + 2S_2,„. {d+i)+ E^2.n2 {d + 2) + 
E2.n2 (d) + 2E2,n2 {d - I) + E^.n^ [d - 2) ^ 8En, . 
By recursion, we can extend this result to any positive index ni 

J2[ )E^n„n2{d + k) + J2[ U;„,,n.(d-fc)-2"i+ii?„,.(A2) 
fc=0 ^""1^ k=0 ^""1/ 

In particular, for n2 = we recover the result ([7]): E-m + £'ni = 2, by using 
the relation J2lLa {nj = 2"^. Equation (gU can be inverted in order to obtain 
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a relation between -E_„i,„2('^) s-'^d functions of positive indices. We use a discrete 
Fourier transform 



to solve the previous Green's function (jA2p 



where we used the Dirac relation 

oo oo 

E ^''"'''= E '^(^-^')- (A3) 



d' — ~oc d' — — oQ 

Inverting this relation, we obtain 



(A4) 



The first integral over the variable z gives simply a Kronecker function (^d-d'-ni.Oj 
whereas, for the second one, we the Dirac sum selects the value z = d' = Q oi z = d' ~ 1 
when performing the integration over z. We obtain the simple result 



/ dz V , , ^ / dz V,5(z-d')7 ^^ = 2- 



-rii 



Therefore, we obtain the symmetry relations between the negative and positive indices 
Til (and 712). Relation ([33| is directly connected to the first two of the lemma, by 
considering the case where the interval lengths are both negative. q.e.d. 

Appendix B. Proof of lemma 3.2 

We consider the case d = in the discrete equation ([Tf]) . Since i?„i^„2(0) = -E'ni+n2J 
we have 

Using ((27)) the first term gives 

C^t^'rai, 712(0) — r 2i?„j+„2 + l — 8i?„j+„2 + -£'ni + l,n2 + -^ni,n2 + l(~l) 

+ 2 En 

1 — 1,712 (l)+i?„,.„2-l(l) , 

whereas the second term gives 

dtEni+n2 — 2r {Eni+n2-l ~ 2£'„j+„2 + ^'ni+n2 + l) ■ 

By comparing the last two expressions, we obtain 

Eni+l,n2{~^) + -E'ni,n2 + l(~l) + -^ni-l,n2(l) + Eni,n.2-l{^) = 4i?„i+„2- 
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By deriving successively the last expression with respect to time, we obtain a general 
Green's equation for the terms involving a negative and a positive distance between 
intervals: 

( ^ Eni+k,n2+d-k{-d) +^ ^ j i?ni -fe,ri2 -d+fe (d) 2''+^£'„j+„2 . 
fc=0 ^ ^ k=0 ^ 

Notice that the sum of all three indices is always equal to ni +77.2. As before, we 
introduce the general Fourier transform 

Enun2{d) ^ / dxdy E{x,y,d)e'^ 
Jo 

00 

E{x,y,d)= Sni.«.(d)e-''"("^^+"^^) 
to obtain the relation 

2d+i 



y) 

rii ,n2 — — 00 



Eix,y,-^d)+E{x,y,d)[ ^ ) = E{x,y,0) 

ier inverse then reads explicitly 
Eni,n2i-d) = y2 -B„'^^„j^(0) / dxdy 

"l,"2 

- Y] Eri'^^n'^id) I dxAy 



The Fourier inverse then reads explicitly 

2d+i 

exp(2i7ra;(?T,i — n'^) + 2i7ry(n2 — n'2)) 

"""" (g2i7rx _|_ g2i7ra)d 

1 / g-2i7rx _|_ g-2i7rj/ ^ d 



g2i7ra: ^ g2i7ry 



g2i7ra;(ni — )+2i7ry(n2 — 



The second integral over x and y gives simply a double Kronecker function 
oi5n2-n' -d,Oj whercas the first integral can be transformed as follows 

V S„;,,4(0) / dxdy 



nj,Ti2 



2d+i 



g2i7r2:(ni — )+2i7ry (n2 —n^) 



^2i7ra:(ni — )+2i7ry (n2 — ?^2) 



"'1 ° 

= V / dxdy 2i^x(n,-n,)+^i^yn2 2i7r«2 (x-y) _ 

n'l ° "2 

The last sum over can be performed using equality (IA3p . which selects x — y over 
the interval of integration over j/. Then, the integration over variable x gives directly 
a Kronecker function 6ni+n2-d-n[,o- Finally we obtain the relation between the index 
—d and d. q.e.d. 

Appendix C. Decomposition of the two-interval probability in the 
three-dimensional space 

We give the detailed calculation for the decomposition of the general solution for the 
two-hole probability. Begin by separating the regions with z' > from those with 
z' < 

Eeo{x,y,z,t) = \ I I dx'dy'dz' We„{x - x' ,y ~ y' , z - z')Eojo{x' ,y' , z') 
+ Wi, [x ~x',y^ y', z + z')E^^t^{x\ y' , -z')] . (CI) 
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In the last term, we map the domain of negative vahies of z' to the positive values by 
using in the continuum limit Eo{x' , y' , —z') — — i?o(x' — z', y'—z', z')-\-2Eo{x' -\-y' ~ 
z'), where Eq{x) = E{x, 0). It is also useful to perform two translations x' — z' ^ a;', 
y' — z' ^ y' in the first term Eq{x' — z' , y' — z' , z'), and a translation on the variable 
x' —i' x' — y' + z' in the second term Eq{x' + — z') so that 



dx'dy'dz' Wia{x -x',y- y',z~ z')Eo,ta{x' ,y' , z') 



Ei^{x,y,z) 



We,{x-x' 



z',y-y' - z',z + z')Eojg (x', y', z') 

+ 2We, ix-x'+y'-z',y-y',z + z')Eo [x'). (C2) 

In the third term containing Eq{x'), we can now perform the integration on y' and 
z' , since the function W^qO is a Gaussian kernel. When x' is negative, we use also 
Eq{x') = 2 — Eo{—x') in order to perform the mapping onto the positive axis. After 
rearranging the different terms we obtain 

Eig (x, y, z) = erfc(z)erfc(a; + y + z) 

da;' 



■ erfc(z) 



-{x+y+z-x')^ 



-{x+y+z+x')^ 



R2 Jo 



dx'dy'dz' 



X - X ,y - y ,z 



- We„ {x-x' - z', y-y' -z',z + z') E„j„{x', y' , z') 



(C3) 



In the last term, the difference between the two weight functions W^qO '^^'^ t)e 
simplified by noticing that 



Wfo ~ x' ,y - y' -.z - z') - Wia {x - x' - z' ,y - y' - z' ,z + z') 



Wi,ix^x',y-y',z-z'){l 



-izz' 



) ■.= WtoW,y',z'). 



(C4) 



From the last multiple integral (|C3p , we can divide the integration domain over x' and 



y' into 4 parts: Eo{x' ,y' , z'), Eq{-x' ,y' , z'), Eo{x' ,-y' , z'), and Eq{- 



-y',z'), with 



(a;', y' , z') all positive variables. Eo{—x', y', z') can be transformed using the following 
steps and and ([M]) 

Eo{-x', y', z') - 2^o(y') - Eo{x', y' , z' - x') 

- 2^o(y') - 0{z' - x')Eo{x', y', z' - x') 

- e{x' - z') VlE^iy' + z') - Eo{z', -~x' + y' + z' , x' - z')] (C5) 

then 

Eoi-x',y', z') = 2E^{y') - 9{x' - z')2Ea{y' + z') 
+ e{x' - z')e{x' -y' - z')2Eo(z') 
-~e{z' --x')Ea{x',y',z' -x') 

+ e{x' - z')e{-x' + + z')Eq{z', -x' + y' + z'., x' - z') 

- 0(x' - z')B{x' -y'- z')Eo{z', x' - y' - z', y'). (C6) 

All arguments of the functions appearing in the last expression are now positive. An 
analogous analysis is done for Eo{x' , ~y' , z'), which is symmetric by inverting x' and 

y' 
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Eo{x', -y', z') = 2Ea{x') - 0{y' - z')2Eo{x' + z') 

+ e{y' - z')e{-x' + y' ~ z')2Eo{z') 
-9{z'-y')Eo{x',y',z'-y') 

+ e{y' - z')e{x' - y' + z')Eo{x' - y' + z' , z' , y' - z') 

- 0{y' - z')ei-x' + y'- z')Eoi-x' + y' - z' , z' , x'). (C7) 

Finally, using (|33|) and (p4| the last expression Eo{—x', ~y', z') is transformed into 
i?o(-.T', -y', z') = 4 - 2Eo{x') - 2Eo{y') + Eo{x' , y', z' - x' ^ y') 

= 4 - 2Eo{x') - 2Eo{y') + 6{z' - x' - y')Eo{x', y', z' - x' - y') 

+ d(x' + y' - z') \2Eq{z') - Eo{z' - y' , z' - x' , x' + y' - z') . 

In the last term, z' — y' and z' — x' can be either positive or negative, which leads us 
to analyse the other possibilities 

Eoi-x', -y', z') - 4 - 2Eoix') - 2Eo{y') + e{x' + y' - z')2Eo{z') (C8) 



+ e{z' - X 

- e{x' + y 

+ 0{y' - z 

+ - y' 

+ 0{y' - 



-y')Eo{x\y\z' -x' ~y') 

- z') \e{z' - y')0{z' - x')Eq{z' - y', z' - .t', x' + y' - z') 

)9{z' — x'){—Eo{y' — z', z' — x , x) + 2Eq{z' — x')} 
)e{x' - z'){-Eo{z' - y\ x' - z', y') + 2Eo{z' - y')} 

)e{x' - z'){A - 2E(i{y' ~ z') ~ 2Eo{x' ~ z') + EQ{y' - z', x' - z', z')} 



Here, we have (i) terms without any empty-interval probability, (ii) terms which refer 
to the single-interval probabilities only and (iii) terms which contain also the two- 
interval probabilities. The first kind of terms is needed for an initially fully occupied 
lattice (see section 4). 



Appendix D. Simplification of some integrals for the initial one-interval 
contribution 

We show how to simplify the integrals Ii,. . . ,Iq of section 4.2. In order to arrive 
at the previous partial result, it is useful to note that a function d{x' — z') (in the 
case when it involves an integral over the one interval distribution Eq which does not 
depend on x') with x' and z' positive can be simplified and replaced by unity if we 
perform the translation x" — x' ~ z' . The integration over x" extends then from 
— z' < to oo, however the function 9 restricts the integration to the interval x" > 0. 
Therefore the limits of the integration are unchanged under this procedure, and the 
function 9 can be set to unity. In the opposite case, where we have a contribution such 
as Eq{x')9(x' — z') for example, a translation on x' will change the argument of Eq, 
which does not conserve the form ([39]) . Instead, we write 9{x' — z') = 1 — 9{z' — x') 
and apply the previous translation to the positive variable z' = z" -I- x' . This changes 
the arguments of the Wiq{---) functions only and not Eq. Using ((37)) . the different 
groups of terms Ii, . . . , I4 can be simplified 

h = mo -v. - w^o W - - 

= (x\ -y\ z') + (x', -y', -z') 
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h = mo i-y', «') - mo i-y' -z',x'-z', z') (Di) 

= m, {-y\ x', z') + mo i-y', x', -z') 

h = mo (-^' - -y' - - mo i-x', -y', z') 

h = Wio i-y' - - - mo i-y', -x', z') 
= - moi-y', -z') - yvto{-y', z'). 

These relations allow us to extend the integration over z' from — oo to +00 since they 
combine each time two terms depending on z' and —z'. Since the variable z' appears 
only in the gaussian weights, the integration over z' gives new gaussian exponentials. 
Then integration over y' gives erf or erfc functions, which are combined with Eq{x'). 
Terms and 1% can be combined together. The first term of 1^, can be transformed 
as 

/ dx'Ay'dz' Eo{x')0{-x' -y' + z')m,{y' ,-z' ,x') 
= ( dx'dy'dz' Eo{x')e{-x' + y' -z')m,{z',-y',x') 

POO POO pO 

= dx'Eoix') dy' dz'6{-x' + y' + z')m,{-z',-y',x'). 

Jo Jo J -00 

The second term gives 

/ dx'dy'dz' Eo{x')e{-x' -y' + z')m,{-z',y',x') = 
Jm.\ 

POO pO poo 

= dx'Eoix') dy' dz'0{-x'+y' + z')m„{-z',-y',x'). 
Jo J -00 Jo 

Finally, consider /g, which is equal to 

pOO POO POO 

/ dx'Eoix') dy' dz'e{-x' + y'-tz')moi-z',-y',x') 
Jo Jo Jo 

and can be combined with the previous two terms. Noticing that 0{—x' + y' + z') =0 

in the domain where y' and z' are both negative, we can write the integral over 

and Iq as follows 

dx'dy'dz' Eo,e,ix')(^h + h) 

= 2V2 r^Eo,,M) ir '^^0{-x' + y' + z')m,{-z',-y',x') 
Jo v'^ JJ-00 

= 2^2 r^Eo,e,ix') [T ^y^e{-x'-^z')m,{y'-z',-y',x'). 
Then, it is possible to perform the integration over y' and z' successively. 

Appendix E. Simplification of some integrals for the initial two-interval 
contribution 

In this appendix, we show how to simplify the expression of the contribution which 
depends on the initial two-interval probability. 
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There are therefore 12 terms involving the weights W^q(- • •)■ After performing different 
translation operations, we arrive at the expression 

E^fix,y,z) \ I Ax'dy'Az' Eo,i,{x' ,y' , z') x 
[m„ ix',y', z') - [^y' - z', ^x' - z\ z') 

+y^lo (-a;', -y', x' +y' + z') - Weg {~x' - z' , -y' - z , x' + y' + z') 
+m,{^x' -y' - z', -z', x' + z') - m,{~x', y', x' + z') 
+W£(,(-2', -x' - y' - z' ,y' + z) - yV(?o(a;', + z') 

+Wt,{-x' - z\ y' + z', x') - i-x' -y' - z\ z\ x') 

+yVtoW + z\ -y' - z',y') - m,iz', -x' - y' - z^y')}. (El) 
We can pair the terms which appear on each line of the previous relation noticing that 
m„ {u~w,v^w,w^{u + v)) ^ {u, v,w~{u + v))e-'^'"^''+y+^l (E2) 
This yields the following relations 





y' 


— z 


1 

~X 
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— z ,z 


)^mo{x\y',z')e-<^ 


+y'+z')ix+y+z) 
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-y ,x +y + z )e 


-42 {x+y+z) 




X 


f 

- y 


1 

— Z 


t 


+ z')=m,i~x\y',x' 
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-x' - 


y'- 


z',y') 


^m,ix' + z',-y'-z' 


y, \^^-Ax' (x+y+z) ^ 





Moreover, we have 

[x', y', z') = Wi, {x -x',y- y', z - z'){l - e^^-'^), 
m„{~^', ~y',x' + y' + z') = W,„(x -x',y^ y',z- ^Oe-^^'^^^^'^il - e-4(-'+!^'+-')-], 
m,X-^',y',x' + z') - W,,{^ -x',y~- y', z ~ 2')e-*^'"[l - e'^^^'+A'], 
W,„(x', ~y',y' + z') = W,„(x -x',y- y',z~ z')e'^y'y[l - e-^^'"' +^'1 
W,, i^x' - z', y' + z', x') =Wi,{x~x',y~y',z- ,>y-^- ix^z) _ ^-4.'.) 

{x' + z', -y' - z', y') = Wi„ {x ~x',y~- y' , z - ^')e- V«-4.'(j/+.) (i _ ^-iy'zy 

Appendix F. Two-point correlation function 

In this appendix, we show how the terms depending on initial conditions are 
corrections to the leading behaviour in the long-time limit. 

We first consider the contribution C^^\z,t) which can be evaluated by deriving the 
kernel in (gOl) 

C^'\z,t) = r E{x\ 0) (erf(z/4)|2 - 4(x' + zf / iAe-^'^'+')' 

- erf(z/4){2 - \{x - z)V^o'}e"("'"")'/^''' (Fl) 
+ 4 r _ ^y^x'-zfl,,,- ^ (2a;' + z)e-^-' ^^^^ I'^^^Xe'^' I'A 
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dx' 

Jo VT^f-o 

In order to obtain the term C^^'^{z,t), we use the important property 
previously seen that ^2(2;', y', 2', 0, 0, z) = dxK2{x' ,y' , z' , x,y, z)\^^y^Q = 
dyK2{x', y', z', X, y, z)]^^^^^ = 0, which can be evahiated from (|Tf| . so that 



C(^\z,t) 



dx' dy' dz' Eo{x', y', z')W{-x\ -y', z - z') x 



- (x' + z')iy' + z')(l + e-4(-'+y')-/^o^y 

'-z'{x' + y' + z') [e-^-'^l'-' + e-V^/^o^^ | /^^4 



da; 







(F2) 



V2 



da;' dy' dz' £0(3;', y', z')W(-x', -y', z - z')L2(x', y', z, z) 



dx ■ 



Ax 



Therefore, we can sum up the different contributions to express the exact correlation 
function. In order to study the effects of the initial conditions in the long-time 
behaviour on the scaling form, we can isolate the dominant contribution for each of 
the three terms C'^"-', C*^^\ C^^\ depending on the properties of the initial distribution 
function for the two intervals and in the case where z/£o is small. For C^^\ we obtain 
easily 



- -^foiz/io), with fo{z/£o) -- 



io 



(F3) 



For C^^\ when £0 is large, we can perform a Taylor expansion of Li(x' , z) around 
x' = 0, and since Li(0, z) = we obtain 



CWiz,t)^il^ rdxEix,0)x 

«0 "'0 



(F4) 



where the scaling function /i is given by the expression 

/i(f ) - ze-^/^«^ \^(^~^) e-'"^' + V^f ( 2^ - 3 ) erf ( f 



z 

lo'»z 16 



4 V 4' 



lb / z" \ 
— 1-6^ 



(F5) 



The quantity Ei :— J^dx E{x,0)x is related to the second moment of the interval 
distribution P{x, 0) (jlOp after performing an integration by parts 



El 



dx E{x, 0) X 



2 7o 



da; P(a;,0)a;2 



(F6) 



Concerning the dominant behaviour of C'^-*, we notice that the gaussian weight 
W{—x', —y', z - z') is peaked around the value (a;', y', z') = (0, 0, z). We can therefore 



C<^'\z,t)c^ I dx'dy%(x',y',z) / .^^L_I_y^(^^' ^ ^y' ^ , ^ ,')L,{x' ,y' , z' , z) 
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begin to perform a Taylor expansfon of Eq{x' , z') around the z' = z and integrate 
the variable z' on the real axis if z is far enough from the value zero, which also is 
satisfied if z/£o small. 

V^dz' 

v ' 

dx'dy'dM^\y',z) f -^^^^{z'-z)W{-x',-y',z~z')L2{x',y',z',z) 
dx^^E{x,0) e-" (F7) 



For £q large, the function Mq can be expanded as a series in x' and y' to obtain the 
dominant term 

The next function Mi in the development can be expanded as 

Mi(x',y',z)^^/«(f) , /«(«) = -«/f (z.). 

iQ <-0 

The last term in (jF7p can be approximated by 

^ Eix,0) e-^/^o^y ^ 

where Ei is the average quantity of the single interval distribution, given in (jF6p . 

Hence, C*-^-* can be expanded as a series of the inverse diffusion length involving 
integrals of the two-interval distribution and related moments, as long as these 
integrals are not diverging: 

C('){z,t)c,l^^]^ f dx'dy'x'y' Eoix',y',z) 



I Ax'dy'x'y' dMx\y\z)-^El. (F8) 

The expansion (|F8[) has the advantage that if the initial condition Eo{x,y,z) does 
not depend on the distance z between the two intervals x and y, then the second 
term of the development involving O^EqIx, y, z) vanishes. However, from the previous 
results (jF3p . (jF4p and (jF8l) . we obtain a hierarchy in the inverse powers l/£o, where 
|C(°'(z,i)| > |C(^)(2;,t)| > |C(^)(z,t)|, with a dominant contribution of order 1/4^, 
1/io and l/(o^i successively. From the previous asymptotic results, it is also easy 
to check that in every case, the sign of the correlation contributions are respectively 
C(°^ < 0, C^^^ > and C*^^^ < 0. In particular, if z < Iq, the dominant contribution 
to C2 comes from the last term of (jF8|) . since the scaling functions /2°'' and f2^'' are 



rapidly decreasing functions of z/£o- Hence 



Ci^H^,t)'"^'^^El (F9) 
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Appendix G. Some integral identities 
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We list some identities involving the kernel W and the error function erf, which are 
used in the main text. 

dx'dy'dz' 



Jo 



^ We„ {x\ y', z') = - erf (x) - erf (x + z) 

o TT 2 L 

dr'dz' - 1 r 

'''' """^ We, {x\ y', z') = - [erf (y) - erf (j/ + z) 



Jo \pK i_oo 

( dr'd7/' Z""" dy' 
4^/2 / / / ^ (x', y', = erf (x) - erf (x + z) 

J Jo ^ -/-oo 

+erf(y) — erf(j/ + 2;) + erf(x)erf (?/) — erf (x + 2)erf(j/ + 2;) 



4V2 ^ 



2\/2 



dy' /• /-^^ dx'dz' 



(a; + + y' - 2;', t/ - t/', a + z) 



= erfc(2;)erfc(x + y + 2) 
dy' ( dx'dz' 



Wi,{-x' -y' + z',y',-z') 



-erf(2;)erfc(x + y + z) 



(Gl) 

(G2) 
(G3) 

(G4) 

(G5) 
(G6) 
(G7) 
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